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The Kagome Heisenberg antiferromagnet is mapped onto an effective Hamiltonian on the star 
superlattice by Contractor Renormalization. Comparison of ground state energies on large lattices 
to Density Matrix Renormalization Group justifies truncation of effective interactions at range 3. 
Within our accuracy, magnetic and translational symmetries are not broken (i.e. a spin liquid 
ground state). However, we discover doublet spectral degeneracies which signal the onset of p6 - 
chirality symmetry breaking. This is understood by simple mean field analysis. Experimentally, the 
p6 chiral order parameter should split the optical phonon degeneracy near the zone center. Addition 
of weak next to nearest neighbor coupling is discussed a . 

PACS numbers: 75.10.Jm, 75.40.Mg 



The antiferromagnetic Heisenberg model on the 
Kagome lattice 

H = J^S Z -S„ J>0, s= 1 -, (1) 

(ij) 

is a much studied paradigm for frustrated quantum 
magnetism. In the classical approximation S —> oo, 
this model exhibits macroscopic ground state degen- 
eracy which encumbers semiclassical approximations. 
There is evidence, both numerical and experimental (in 
ZnCu3(OH)eCl2 [1]), that quantum fluctuations lead to 
a paramagnetic "spin liquid" ground state [2]. 

Exact diagonalization studies (ED) [3], and ED in 
the variational dimer singlets subspace [4], have not ap- 
proached the thermodynamic limit, due to severe com- 
puter memory limitations. Many methods have proposed 
paramagnetic ground states, including lattice symmetry 
breaking "valence bonds crystals" [5-8], algebraic spin 
liquids [9] and a time reversal symmetry breaking, chiral 
spin liquid [10]. 

To date, the lowest energy on long cylinders has 
been found by Density Matrix Renormalization Group 
DMRG [11, 12]. The DMRG ground state is a trans- 
lationally invariant singlet, with apparently no broken 
translational or rotational symmetries. This state is con- 
sistent with a resonating valence bonds (RVB) state [13] 
with a spin gap As=i = 0.13, (henceforth we express en- 
ergies in units of J), and Z2 topological order [12, 14]. It 
is still unclear however, what are the low-energy singlet 
excitations of this state [11, 15], and whether or not any 
other symmetry of % may be broken in the infinite two 
dimensional limit. 



a This work was supported by the U. S. DOE, Contract No. DE- 
AC02-76SF00515. 




FIG. 1. The CORE-blocking scheme on the Kagome lattice. 
R and L denote the two pinwheel ground states of the 12 site 
stars, and the arrows (pseudospins) denote the symmetrized 
Ising basis which spans the reduced Hilbert space of }j CORE ? 
see Eq. (2). 

This paper reports a surprising result: The thermody- 
namic ground state appears to break reflection symme- 
tries, and to possess two dimensional p6 chirality, (- not 
to be confused with "spin chirality" which also breaks 
time reversal symmetry [10]). Our conclusion is obtained 
by Contractor Renormalization (CORE) [16] with 12-site 
stars blocking, see Fig. 1. the stars scheme is found to 
reach sufficient accuracy with range-3 interactions. This 
is evidenced by comparing ground state energies of the 
effective Hamiltonian H COREs , to high precision DMRG 
on large lattices. The small modulation of bond energies 
is consistent, within our accuracy, with a translationally 
invariant singlet state as deduced by DMRG [11, 12]. 

■ffCORE 3 - g diagonaiized Q n up to 27 stars (effectively 
324 Kagome sites). The spectra exhibit doublet degen- 
eracies between states with opposite parity under reflec- 
tion [17]. These signal a hitherto unexpected sponta- 
neous symmetry breaking in the thermodynamic limit 
into a chiral ground state. This chirality is understood as 
the effect of three-star interactions. Classical mean field 
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theory on the effective hamiltonian explains this symme- 
try breaking. A two-dimer chirality order parameter is 
defined on the microscopic Kagome model. We propose 
an experimental signature of this broken symmetry in 
the phonon spectrum: a splitting of symmetry-protected 
degeneracy between two zone center optical modes [18]. 

Finally, we add weak ferromagnetic next nearest neigh- 
bor interactions J2, and find that it eliminates the chi- 
rality at J 2 ~ —0.1. 

CORE procedure. Previous CORE calculations for the 
Kagome model [19, 20] started with up-triangles block- 
ing, and did not reach sufficient convergence at range 3. 
Here we use much larger and more symmetric blocks of 12 
site (Magen David) stars which form a triangular super- 
lattice. In each star, we retain just the two degenerate 
singlet ground states \Li) and \Ri), depicted in Fig. 1 
which form a pseudospin-1/2 basis: 



ti> = 
\u) 



^2+1/16 

-J 

V/2-1/16 



\Ri) 



\Ri) 



\Li)), 
\Li))- 



(2) 



Note that the two states are invariant, and have op- 
posite parity under all Dq reflections. 

The CORE effective Hamiltonian on a superlattice of 
size N s stars is defined by the cluster expansion, 



N s 



ff COflE = E^ (1) +E + E d 

i=l i\i<2 ii,...ijv 

e^^-E E 4 

m<n f3(m)£a(n) 



Am) 



(3) 



Here f3(m) is a connected subcluster of size m in a clus- 
ter a(n) of size n stars, and is defined to be an in- 
teraction of range-n. The operators are constructed 
by ED of Eq. (1) on a Kagome cluster a: 



2 n 

E e "i*")^"i 



(4) 



Here (e",\I/") are the exact 2 n lowest singlet energies 
and wavefunctions. The states are an orthogonal 

basis constructed by sequential projections of = 
1, 2, . . . 2 n onto the pseudospin states. After projection, 
the states are orthogonalized sequentially by using the 
Gram- Schmidt procedure. 

If interactions of all ranges n < N s are included, then 
H CORE has the identical low energy singlet spectrum as 
Eq. (1) on the equivalent Kagome lattice. However, ED 
cost to compute grows exponentially with n. Thus, 
the success of CORE depends on the ability to truncate 
the cluster expansion at feasible n while maintaining suf- 
ficient accuracy in the truncated Hamiltonian. 

The error in the ground state energy 8Eq OREji can 
be computed by comparison to high-precision DMRG 
on large lattices with m > n stars. This error should 
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J x 


J y 


J z 


ED 


-6.26391 


0.13818 


0.00713 


-0.00105 


-0.00045 


PT 


-5.268 


0.046 





-0.00025 


-0.00175 



TABLE I. Parameters of CORE range-2 Hamiltonian, by 
Exact Diagonalization, and second order Perturbation The- 
ory [22]. 



be much smaller than the important interactions in 

jjCORE n 

Lattice translations. Our choice of stars for the reduced 
Hilbert space nominally breaks lattice translational sym- 
metry as seen in Fig. 1. The microscopic spin correlations 
are computed by functional differentiation of the CORE 
ground state energy with respect to source terms [21]. In 
principle one must compute the effective interactions to 
all ranges to restore full translational symmetry. Nev- 
ertheless, symmetry breaking artifacts decrease with the 
truncation range n. We can therefore identify any spon- 
taneous translational symmetry breaking which signifi- 
cantly exceeds the truncation error. 

CORE range 2. We start with the lowest-order trun- 
cation at range 2. The general form of the two-star in- 
teractions allowed by lattice reflection symmetries is 



H CORE *=Nc + hJ2^ 



+ E^ 

(ij) 



a —.ex -.ex 

(J A <J A , 



(5) 



where i labels sites, (ij) nearest neighbor bonds on the 
triangular lattice. cr a ,a = x,y,z are Pauli matrices. 

The parameters derived from the lowest 4 eigenstates 
of 24 spins, are computed by Lanczos algorithm, and 
listed in Table I. It is instructive to compare the ED pa- 
rameters to the second order perturbation theory (PT) in 
the inter-star bonds, as was calculated by Syromyatnikov 
and Maleyev citeSM. Second order PT in the connecting 
bonds is not very accurate when connecting bonds have 
exchanges equal to 1. For example, PT misses the im- 
portant J x interactions. The dominant interaction of 
h core 2 ig the field h = 0.138, which would yield in the 
thermodynamic system a ferromagnetic ground state po- 
larized in the | \) direction. In terms of Kagome spins, 
the ground state would be a product of antisymmetric su- 
perposition of pinwheel states, with local t fluctuations 
generated by the xx, yy terms. 

Within CORE2, the connecting bonds energy is 
Einter = —0.21283, versus the intra-star bonds at 
E i n t r(l = —0.2225. Interestingly, the modulation is al- 
ready diminished from 100%-^4.3% with range-2 inter- 
actions. 

How accurate is H CORE2 ? Unfortunately, not 
enough. The exact ground state energy/site of H for 
36 sites is E ED =-0.41276 while the CORE 2 energy/site 
for three stars is E° ORE2 =-0.4277. The error of -0.0149 
is quite large relative to the important terms in H CORE2 . 
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J 2 


J 2 — (Kagome) 


J 2 = +0.1 


J 2 = -0.1 


CO 


-5.24629 


-5.17068 


-5.48631 


h 


-0.069224 


0.059323 


-0.362797 


J X 


-0.009028 


-0.015421 


0.001123 


Jy 


-0.011879 


0.001832 


-0.017699 


Jz 


0.021056 


0.003686 


0.020141 


J zxx 


-0.027920 


-0.019649 


-0.009524 


Jzyy 


0.004550 


-0.004749 


0.002394 


J zzz 


0.000660 


-0.001410 


0.010495 



TABLE II. Interaction parameters of CORE range 3, with 
three values of Ji. 



Hence we cannot trust the truncation, and longer range 
interactions are needed. 

CORE range 3. To obtain H COREs we compute the 
interactions on the three star triangular cluster [23]. 
This required ED of Eq. (1) of 36 spins with open bound- 
ary conditions (OBC). For verification, we ran both a 
standard Lanczos routine on a supercomputer, and the 
memory-economical Lanczos- SVD routine [24] on a desk- 
top computer. Adding contributions from ranges 1- 3 we 
obtain the following effective hamiltonian 



H CORE 3 = Ncq 



haf 



^2 Ja(J i 

(ij),a 



(6) 



where (ijk)^ label nearest neighbor triangles on the tri- 
angular lattice. We do not list terms that cancel in the 
superlattice summation with Periodic Boundary Condi- 
tions (PBC). 

The H CORE z truncation error is estimated by compar- 
ing its ground state energy on clusters of sizes TV > 3 to 
high precision DMRG [25, 26]. From Table III we find 
that the truncation error is very satisfactory < 0.004 per 
site. If we extrapolated CORE3 ground state energy to 
the thermodynamic limit, we get —0.447 which underes- 
timates the extrapolated DMRG result -0.439 [12] by at 
most —0.008 per site. We see in Table II for comparison, 
that the values of the dominant terms /i, J zxx and J z 
(multiplied by number of bonds/site), are significantly 
larger than this error. Thus, we believe that CORE trun- 
cation at range 3 is sufficiently accurate to obtain the 
correct thermodynamic phase. 

p6-Chirality. The ED spectrum of H COREs is evalu- 
ated on lattices of up to N s = 27 stars (324 Kagome 
sites), with PBC. The most striking feature on lattices 
larger than N s = 9, is the emergence of ground state 
degeneracy of two singlets with opposite parity under re- 
flections. In the pseudospin representation, even (odd) 
parity states include only an even (odd) number of stars 
with antisymmetric | \) states. These degeneracies signal 
a spontaneous reflection symmetry breaking p6m — >• p6 



number of stars 


e core 3 


DMRG 


Error 


2x2 


-0.418452 


-0.417213 


-0.001239 


2x3 


-0.423953 


-0.422336 


- 0.001617 


3x4 


-0.431150 


-0.428046 


-0.003104 


3x5 


-0.432688 


-0.429191 


-0.003497 



TABLE III. Ground state energies per site of H CORE3 , and 
comparison to DMRG on equivalent Kagome clusters (with 
OBC). 




FIG. 2. a) The mean field ground state of H CORE3 exhibiting 
(a x ) > order, which corresponds to two dimensional chiral- 
ity. b) A typical singlets configuration in the corresponding 
ground state of the Kagome lattice. Notice that there is no 
translational order, but that there are more pinwheel con- 
figurations of \R) than \L). c) The two dimer chiral order 
parameter defined in Eq.8. 



in the thermodynamic limit. 

A Mean Field (MF) energy of H CORE3 in spin- 1/2 
coherent states |f^) is 

E MF =7Vc + /i^cos6> i + •W" n ?> 

+ Yl JzaaCOS^n?^, (7) 

{ijk),a 

where = (sin 9i cos 0^, sin 0i sin cos 0j). In Table II 
we see that for J 2 = 0, the dominant couplings are the 



J2 


M MF 


Mf D 


0.0 


0.2647 


0.2257 


+0.1 


0.1390 


0.1476 


-0.1 


0.5 


0.4999 



TABLE IV. Effect of next nearest neighbor couplings on the 
ground state ^-polarization. 



field h, and the J z and J zxx exchanges. The last coupling 
is responsible for the chiral symmetry breaking, as it pulls 
the spins in the ±£ direction. 

Minimizing E MF ', we find a ferromagnetic state de- 
picted in Fig. 2(a). The z-polarization M^ F = ^cos(9 
is compared to ED in Table IV. For J 2 = 0, we find that 
the chirality order is substantial with ^ sin# = 0.424 by 
MF, and \{a x ) = 0.397 by ED. 

In Fig. 2(b) we depict a typical dimer configuration 
which contributes to the p6 - chiral RVB state. One can 
see the predominance of R pinwheel chirality over L. The 
most local order parameter for this chirality is the two 
dimer correlation depicted in Fig. 2(c), 

Ci = ^2 (^ d ^V(d) — £d<^V(d)) ? (8) 

d 

where the dimer singlet projectors are 

S d = l/4-S dl .S da , (9) 

and ?7 r (d) (r/(d)) is the bond emanating from i at angle 
7r/3 (27r/3) relative to the dimer bond opposing i. The 
two terms in C measure parts of pinwheels of opposite 
chirality. 

Translational symmetry. At range 3, the energy of 
internal triangles E/± inter — —0.686, and connecting tri- 
angles is E/±. ntra — —0.665 (depicted by solid and dashed 
lines respectively in Fig. 1). This relative modulation of 
about 3.0% lies within the truncation error. Thus we can 
affirm that CORE 3 ground state is consistent with trans- 
lational invariance in agreement with DMRG [11, 12]. 

Singlet Excitations. In the 27 star lattice, the low- 
est singlet excitation above the two degenerate ground 
states is AEs=o = 0.28, which has a non zero wavevec- 
tor. This excitation gap does not vary much with lat- 
tice size. Within the pseudospin Hamiltonian, it can be 
understood as a local spin flip from the ferromagnetic 
ground state. We note that the singlet gap is slightly 
higher than two S = 1 magnons at energies Es=i = 0.13. 
This conclusion differs from that obtained by ED on 36 
site PBC, which found a large number of singlets be- 
low the spin gap [27]. Since our effective Hamiltonian 
describes excitations on much larger lattices, we are in- 
clined to associate these low singlets with the smaller 
PBC lattice geometry. 

Experimentally, fluctuating two-dimer correlations are 
tricky to observe directly. Fortunately, real compounds 
have sizable magneto-elastic coupling between the ions 
and the dimer singlets. While, on average, dimer density 
and bond lengths are uniform in the RVB state, dimer 
density fluctuations, Spd governed by the characteristic 
singlet energy scale, are linearly coupled to the ionic dis- 
placements. In Fig. 3, the effect of a temporary excess 
of dimers on a triangle is shown. In the chiral phase, 
imbalance between the left and right bonds emanating 
out of the triangle produces a chiral force on the ions as 
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FIG. 3. Kagome phonon spectra in p6m phase, and p6 
- chiral phase, calculated within a nearest neighbor spring 
constant model given in Ref.[28]. In the top right, the dimer 
chiral correlations induce a linear coupling between excess 
dimer density S pa and chiral ionic displacements, as depicted 
by the arrows. This adds a chiral term to the dynamical 
matrix which splits the degeneracy of the optical phonons at 
the zone center. 

depicted by the arrows. Integrating out the dimer den- 
sity fluctuations results in a chiral perturbation to the 
phonon dynamical matrix [17]. By symmetry [18], the 
degeneracy between two optical modes is removed at the 
zone center, as shown in Fig. 3. 

Finite J 2- We have added next nearest neighbor in- 
teractions with coupling J 2 to Eq. (1), and calculated 
the parameters of H COREs , as shown in Table II. For 
J 2 = 0.1, we find the same doublet degeneracies, and 
chirality, as for the pure model J2 = [14]. In contrast, 
for a weak negative J 2 = —0.1 the spectrum changes 
dramatically: The doublets are removed, and the ground 
state is fully polarized in the t direction. The precise 
nature of this phase has not been yet explored. Interest- 
ingly, we notice that in proximity to the parameters of 
Table II, one finds the Ising antiferromagnet in field. Its 
ground state contains ferromagnetic hexagons with re- 
versed spins in their center. It represents the Hexagonal 
Valence Bond Solid state, previously shown to have low 
variational energies [5], and proposed for J2 — —0.1 [29]. 

Summary. Using CORE we arrived at an effective 
Hamiltonian, whose accuracy was determined to be suf- 
ficiently high so as to trust its predictions for the ther- 
modynamic limit. Its ground state is consistent with a 
translationally invariant RVB phase, but with broken p6 
chiral symmetry. A 2-dimer chiral order parameter is 
defined, which may be numerically explored on large lat- 
ices. Experimentally, it may be detected by splitting of 
optical phonon degeneracy. 

Acknowledgements. We thank Andreas Lauchli, Daniel 
Podolsky and Didier Poilblanc, for useful discussions. 
A A and SC acknowledge the hospitality of the Aspen 



5 



Center for Physics and KITP at Santa Barbara, sup- 
ported by the NSF Grants PHY- 1066293 and PHY- 
1125915. Funding from U. S. DOE, Contract No. DE- 
AC02-76SF00515, Israel Science Foundation and U.S. - 
Israel Binational Science Foundation are acknowledged. 
Numerical simulations were performed at CALMIP. 



[1] Oren Ofer et al, (unpublished) arXiv:cond-mat/06 10540; 
J. S. Helton et a/., Phys. Rev. Lett. 98, 107204 (2007). 
A. Olariu, et al, Phys. Rev. Lett. 100, 087202 (2008); 
K. Matan et a/., Phys. Rev. B 83, 214406 (2011). 

[2] L. Balents, Nature (London) 464, 199 (2010). 

[3] C. Zeng and V. Elser, Phys. Rev. B 42, 8436 (1990); 
P. Lecheminant, B. Bernu, C. Lhuillier, L. Pierre, and 
P. Sindzingre, Phys. Rev. B 56, 2521 (1997); C. Waldt- 
mann, H.-U. Everts, B. Bernu, C. Lhuillier, P. Sindzin- 
gre, P. Lecheminant, and L. Pierre, Eur. Phys. J. B 2, 
501 (1998); H. Nakano and T. Sakai, Jour. Phys. Soc. 
Japan, 80, 53704 (2011); A. M. Lauchli, J. Sudan, and 
E. S. S0rensen, Phys. Rev. B 83, 212401 (2011). Recently, 
48 spins have been diagonalized (A. M. Lauchli, private 
communication) . 

[4] C. Zeng and V. Elser, Phys. Rev. B 51, 8318 (1995); M. 
Mambrini and F. Mila, Eur. Phys. J. B 17, 651 (2000). 

[5] J. B. Marston and C. Zeng, J. Appl. Phys. 69, 5962 
(1991); P. Nikolic and T. Senthil, Phys. Rev. B 68, 
214415 (2003); R. R. P. Singh and D. A. Huse, Phys. 
Rev. B 77, 144415 (2008). 

[6] G. Evenbly and G. Vidal, Phys. Rev. Lett. 104, 187203 
(2010). 

[7] D. Poilblanc, M. Mambrini, D. Schwandt, Phys. Rev. 

B 81, 180402(R) (2010); D. Poilblanc and G. Misguich, 

Phys. Rev. B 84, 214401 (2011). 
[8] B. Clark, J. Kinder, E. Neuscamman, G. K.-L. Chan, M. 

J. Lawler, arXiv: 1210. 1585, unpublished. 
[9] M. Hermele, Y. Ran, P. A. Lee, and X. G. Wen, Phys. 

Rev. B 77, 224413 (2008); Y. Iqbal, F. Becca, and D. 

Poilblanc, Phys. Rev. B 84, 020407 (2011). 
[10] L. Messio, B. Bernu, and C. Lhuillier, Phys. Rev. Lett. 

108, 207204 (2012). 



[11 

[12; 

[13 
[14 

[is; 

[16 



[17 

[is; 

[19 

[20; 

[21 

[22; 
[23; 

[24; 

[25 
[26' 

[27 

[28; 
[29; 



S. Yan, D. A. Huse, and S. R. White, Science 332, 1173 
(2011). 

S. Depenbrock, I. P. McCulloch, and U. Schollwock, 

Phys. Rev. Lett. 109, 067201 (2012). 

P. Fazekas and P. W. Anderson, Philos. Mag. 30, 423 

(1974). 

H.-C. Jiang, Z. Wang, and L. Balents, arXiv:1205.4289. 
H. C. Jiang, Z.Y. Weng, and D. N. Sheng, Phys. Rev. 
Lett. 101, 117203 (2008). 

C. J. Morningstar, M. Weinstein Phys. Rev. D 54, 4131 
(1996); M. Weinstein, Phys. Rev. B 63, 174421 (2001); 
M. S. Siu and M. Weinstein Phys. Rev. B 77, 155116 

(2008) . 

Supplementary material. 

We are indebted to Daniel Podolsky for suggesting to us 
this idea. 

R. Budnik and A. Auerbach Phys. Rev. Lett. 93, 187205 
(2004). 

S. Capponi, A. Lauchli, and M. Mambrini Phys. Rev. B 
70, 104424 (2004). 

E. Altman and A. Auerbach, Phys. Rev. B 65, 104508 
(2002); E. Berg, E. Altman, and A. Auerbach, Phys. Rev. 
Lett. 90, 147204 (2003). 

A. V. Syromyatnikov and S. V. Maleyev, Phys. Rev. 
B 66, 132408 (2002); A. V. Syromyatnikov and S. V. 
Maleyev, JETP 98, 538 (2004). 

We order the clusters in terms of their radii. Thus three 
star lines are considered longer range than the triangle. 
Three star 120° are not included in the cluster expansion, 
as they get cancelled in embedding range 4 rhombi [21], 
and E. Altman, Ph.D. thesis, Technion, unpublished. 
M. Weinstein, A. Auerbach, and V. R. Chandra, Phys. 
Rev. E 84, 056701 (2011). 

S. R. White, Phys. Rev. Lett. 69, 2863 (1992). 

We use open boundary conditions on these clusters, and 

keep up to to get a discarded weight smaller than 10 -5 . 

Note that these variational DMRG energies give exact 

upper bounds on the ground-state energy [11]. 

P. Sindzingre, C. Lhuillier, Europhys. Lett. 88, 27009 

(2009) . 

L. Zhang, J. Ren, J.-S. Wang and B. Li, J. Phys.: Con- 
dens. Matter 23, 305402 (2011). 

O. Ma and J. B. Marston, Phys. Rev. Lett. 101, 027204 
(2008); Y. Iqbal, F. Becca, and D. Poilblanc, Phys. Rev. 
B 83, 100404(R) (2011). 



Supplement: p6-Chiral Resonating Valence Bonds in the Kagome Antiferomagnet 



Sylvain Capponi 1 , V. Ravi Chandra 2 ' 3 , Assa Auerbach 3 , and Marvin Weinstein 4 
1 Laboratoire de Physique Theorique, Universite de Toulouse and CNRS, UPS (IRSAMC), F-31062, Toulouse, France 

2 Physics Department, Technion, Haifa 32000, Israel 
3 National Institute of Science Education and Research (NISER), IOP Campus, Bhubaneswar, 751005, India and 
4 SLAC National Accelerator Laboratory, Stanford, CA 94025, USA. 
(Dated: October 29, 2012) 

Here we provide details of the exact diagonalization of h core 3 on 27 stars, and the effects of the 
chiral order parameter on the optical phonons of the Kagome ions. 
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EXACT DIAGONALIZATION OF 27 STARS 

In Fig. 1 we plot the exact spectrum of 27 stars with 
periodic boundary conditions. Since this cluster has all 
translations and the full Cq v point group symmetry, we 
can label eigenstates according to their symmetry quan- 
tum numbers. The 6 inequivalent momenta are shown in 
the plot, and for each of them, we can make use of the 
little group, i.e. the reduced symmetry group that leaves 
the momentum invariant: at T point, the little group is 
Cq v , at K point C% v , while A, B, C, and D momenta only 
have a single Z2 reflection. Therefore, we can label few 
low-energy states according to the irreducible represen- 
tations of these groups. 
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FIG. 1. Top: a cluster of 27 sites on the triangular super- 
lattice, and its wavevectors in the Brillouin Zone Spectrum. 
Bottom: The singlet spectrum of H CORE 3 on this cluster 
(equivalent to 324 Kagome spins). Notice the nearly exact 
doublet degeneracies, including the ground state doublet at 
TAi, FA2. These degeneracies signal p6m — »• p6 spontaneous 
symmetry breaking in the thermodynamic limit. 



COUPLING OF THE CHIRAL ORDER 
PARAMETER AND THE PHONON SPECTRUM 

p6m Model 

A simple model for the phonon spectrum of the planar 
Kagome was worked out by L. Zhang et al. in Ref. [1]. 
The two dimensional ionic dispacements are labelled 
by i =1= bottom left, 2= bottom right, and 3= top site. 
The nearest neighbor force constant matrix in the xy 
basis is 



K L 
K T 



(1) 



where we Kl and Kt are longitudinal and transverse 
force constants respectively. We follow Ref. [1] and take 
K T = Kl/4. A rotation is defined: 



U 



cos(7r/3) — sin(7r/3) 
sin(7r/3) cos(7r/3) 



(2) 



with which we express the rotated force constant matrices 

K 12 = K X , K 13 = UK X U\ K 23 = U^K X U (3) 

These matrices define blocks in the 6x6 spring constant 
matrices 



K n = 



K 3 




-K 12 -K 13 
2(K 12 + K 23 ) -K 23 
-K 23 2(K 13 + K 23 ) 


K 2 = I 

-if 13 

-K 12 
K 4 = I 












-K 23 




(4) 



K 2 e^+^ 



The dynamical matrix in q-space is 
K (q) = K + K x e iq * 

+K 3 e i{ ~^ + ^ ) +K 4 e- iq - 



y/3q y 



(5) 
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FIG. 2. Kagome phonon spectra in p6m and p6-Chiral phases, 
calculated within a nearest neighbor spring constant model 
given in Ref. [1]. In the top right, the dimer chiral correlations 
induce a linear coupling between excess dimer density and 
chiral ionic displacements, as depicted by the arrows. This 
adds a chiral term to the dynamical matrix and splits the 
degeneracy of the optical phonons at the zone center. 

Chiral Symmetry Breaking RVB - A magneto-elastic 
coupling |^ exerts a force on the atoms on either side of 
a singlet dimer toward each other. However, the RVB 
phase has a uniform average density of dimers on all 
bonds. There are however dynamical fluctuations of the 
dimer density which can be modelled by a harmonic La- 
grangian 

C d =\f d 2 x\{d 2 t -ul)5p 2 d {x) (6) 
where Ud is of the order of the singlet excitation energies. 



Let us consider a RVB phase with positive chirality 
given by (C) > 0. An instantaneous excess dimer density 
on the up-triangle is shown in Fig. 2. (In any dimer cov- 
ering, some fraction of the triangles have no dimers, and 
thus lower than the average density). Due to right chi- 
rality, the neighboring triangles' dimers directions are bi- 
ased to emanate rightward relative to the opposite edge. 
As a consequence, each corner of our triangle feels a net 
"pull" along the direction of its (on the average) uncom- 
pensated bond. This amounts to an effective coupling 

C d - lon = - 7 |(C> | J d 2 x5p 2 d (x)v r • u(x) (7) 

where u is the ion displacements vector, and v r is given 
as the 6-dimensional clockwise displacements, 

v r = -^( c °s(|),sin(|),-l,0,cos(|),-sin(|) 

By integrating out the gapped dimer density fluctuations, 
a chiral correction to the quadratic dynamical matrix is 
obtained 

D chiral (ci) = D(q) - 72 !^V v r (8) 

This correction lifts the degeneracy of the optical 
phonons at q = as shown in Fig. 2. 



[1] L. Zhang, J. Ren, J.-S. Wang and B. Li, J. Phys.: Condens. 
Matter 23, 305402 (2011). There is a typographical error 
in the labeling of the Ki matrices was corrected here. 



